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Abstract 

We present a study of electroweak production of top and antitop quarks in the s-channel mode 
at the LHC, including next-to-leading order (NLO) quantum chromo dynamics (QCD) corrections 
to the production and decay of the single (anti)top quark. The spin is preserved in production 
and decay by using the narrow width approximation for the (anti)top quark. We show the effect 
of different 0(a s ) contributions on the inclusive cross section and various kinematic distributions 
at parton level after imposing relevant kinematic cuts to select s-channel single top quark events. 
We also discuss several possibilities for measuring the top quark polarization. 
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(a) (b) W 

Figure 1: Representative Feynman diagrams of the three single top quark production modes: t- 
channel (a), associated production (b) and s-channel (c). 

I. INTRODUCTION 

Recent results at the Tevatron pp collider have confirmed the existence of electroweak 
single top quark production l|-|5|. While the Tevatron can be considered a ti factory, 
measurements of single top quark properties are statistics limited. This limitation does 
not exist at the LHC The main mode of top quark production at the LHC is still strong 
interaction top quark pair production, but the number of produced single top quark events 
will be large enough for precision measurements. 

Single top quark events are of considerable importance for probing the Standard Model 
of particle physics (SM). As the top quark decays via the weak interaction before it can 
hadronize, it is possible to measure its polarization. In single top quark events, the top 
quark is coupled to the bottom quark with an amplitude proportional to the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element Va, so that a value for Vt& can be obtained by 
measuring the single top quark production cross section. 

Electroweak single top quark production at the LHC occurs in three different modes (cf. 
Fig. [I]). The process with the largest cross section is the t-channel exchange of a virtual W 
boson (bq — > tq' and bq' — > tq), also referred to as IV-gluon fusion, followed by associated 
production of a top quark and a W boson (bg — > tW~) and the s-channel decay of a virtual 
W boson (qq' ->■ W* -> tb). 

The top quark is the heaviest known elementary particle, which makes it an excellent can- 
didate for new physics searches. The cross sections of the three single top quark production 
modes are sensitive to different kinds of new physics 

Q-Q- The 

single top quark s-channel 
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cross section is especially sensitive to additional bosons, which makes it a very important 
channel in spite of its comparably small cross section. The t-channel provides good means 
to search for flavor changing neutral currents (FCNC) and the associated production cross 
section changes with a modified Wtb or Wqq' coupling. Both t-channel and Wt associated 
production are sensitive to the bottom quark parton distribution function (PDF). A detailed 
knowledge of the properties of single top quark production is necessary, as it is an impor- 
tant background for several Higgs boson production modes. The Wt associated production 
serves as background to Higgs boson searches in the decay channel H — > WW 26| and also 
to processes involving a charged Higgs boson like bg — > tH ± and H ± — > rv 2^] . 

In order to extract the single top quark signal from the large QCD and W + jets back- 
grounds, but also in cases in which single top quark production is a background itself, accu- 
rate theoretical predictions including higher order QCD corrections are needed. The NLQ 
QCD corrections to the single top quark production have been carried out in Refs. 
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37]. 



Furthermore, the complete NLO QCD calculations including both single top qu ark produc- 
tion and decay have been carried out for the Tevatron in several studies 38M43| . To a lesser 
extend, single top quark production is also affected by electroweak corrections and possibly 



virtual supersymmetric effects 



44N46| . In this paper we use the full NLO QCD calculations 



for single top quark production at the LHC in order to present a detailed phenomenological 
analysis, focusing on signal cross sections and kinematical distributions at parton level after 
imposing simple kinematic cuts. 

In contrast to the Tevatron, the LHC is a pp collider and we have to analyze top and 
antitop quark production separately. In this paper, we present distributions for top quark 
production alone where the top-antitop quark differences are small and contrast them with 
the results for the antitop quark where they are not. 

In Sec. [Til we first present the inclusive cross section for s-channel single top quark 
production and discuss its dependence on the center of mass energy of the collider (E c , m _), 
the top quark mass (mt) and renormalization and factorization scales. We also evaluate PDF 
uncertainties. In Sec. we examine the effect of various kinematic cuts. Kinematical 
distributions of final state objects and spin correlations are discussed in Sec. IIVI Our 
conclusion follows in Sec. |V] 
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II. CROSS SECTION (INCLUSIVE RATE) 



In this section, we show the inclusive production rates for s-channel single top quark 
production and discuss their dependence on E c m , m t and factorization and renormalization 
scales. We present numerical results for s-channel single top quark events considering the 
leptonic decay of the W boson from the top quark decay at the LHC (a pp collider). Unless 
otherwise specified, we use the NLO parton distribution function (PDF) set CTEQ6.6M 
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51] , defined in the MS scheme, and the NLO (2-loop) running coupling a s with Ajj^ provided 



by the PDFs. For the CTEQ6.6M PDFs, = 0.326 GeV for four active quark flavors. 

Ad o 

For the numerical evaluation, we choose the following set of SM input parameters: = 
1.16637 x 10~ 5 GeV" 2 , M w = 80.413 GeV, M z = 91.187 GeV, and a s (M z ) = 0.1186. The 
square of the weak gauge coupling is g 2 = Ay/^M^G^. 

In this study we focus on the electron leptonic decay of the W boson from the top quark 
only, but for muon leptons the analysis procedure would be analogue. Including the 0(a s ) 
corrections to W — > qq'. the branching ratio for the decay of the W boson into leptons is 



Br(W l + v) = 0.108 52] . If not otherwise specified, the top quark mass is chosen to be 



m t = 175 GeV 53[, the center of mass energy of the collisions E c _ m _ = 14 TeV and we will 



choose the renormalization scale (hr) as well as the factorization scale (/ip) to be equal to 
m t . In the current section we present inclusive cross sections, which include all W boson 
decay modes. 

In order to calculate NLO QCD differential cross sections we adopt the one-cutoff phase 
space slicing (PSS) method [54|-l56j with a cutoff parameter s min = 5 GeV 2 . 



Inclusive cross section 



As in our previous studies 40|, |41[, we divide the higher-order QCD corrections into 
three separate gauge invariant sets: corrections to the initial particles (INIT), corrections 
to the final state (FINAL), and corrections to the top quark decay (SPEC). The explicit 



diagrams and definitions for the different corrections can be found in Ref. [39j. For E c 
14 TeV and m t = 175 GeV, Table [J shows the inclusive cross sections for top and antitop 
quark production as well as the individual 0(a s ) contributions. The effects of the finite 
widths of top quark and W boson have been included. The total NLO s-channel single top 



4 



Top Antitop 





Cross section 


Fraction of 


Cross section 


Fraction of 




[Pb] 


NLO (%) 


[pb] 


NLO (%) 


Born level 


4.42 


72.55 


2.70 


72.07 


INIT 


1.18 


19.38 


0.73 


19.48 


FINAL 


0.80 


13.16 


0.51 


13.57 


SDEC 


-0.31 


-5.09 


-0.19 


-5.11 


0(a s ) sum 


1.67 


27.44 


1.04 


27.93 


NLO 


6.09 


100 


3.74 


100 



Table I: Inclusive single top quark production cross sections for different subprocesses, for top 
quark production (left) and antitop quark production (right). -E c . m . = 14 TeV and mt = 175 GeV. 



281 ]. but updated values for the electroweak 



quark production cross section agrees with Ref. 
parameters are used. 

As can be seen in Table HJ the LO cross section for antitop quark production (ud — > ib) is 
39% smaller than the cross section for top quark production (ud — > tb). This is due to the 
difference in parton densities of the colliding protons. While in both cases the antiquark is 
from the quark sea of one of the incoming protons, the probability that it collides with an 
up quark from the other proton is higher than the probability for a collision with a down 
quark. 

The 0(a s ) corrections increase the cross section by 38% for both top and antitop quark 
production. The largest contribution to the 0(a s ) corrections comes from the initial state, 
due to the enhancement of collinear physics and the large phase space for additional parton 
radiation. As expected, and as also observed at the Tevatron, the correction to the top 



quark decay is small 



40 



4l|. 



The obtained inclusive cross sections at the LHC for s-channel single top quark events 
(considering the W boson decay branching ratio and using m t = 175 GeV, 172.5 GeV and 
170 GeV) are shown in Table HT1 and Fig. [2] for top and antitop quark production at different 
E c .m.- For the LHC injection energy, the cross sections become very small, for m t = 175 GeV 
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Top 


Antitop 


m t [GeV] 


Ec.m. [TeV] 


Cross section [pb] 


Cross section [pb] 




14 


6.09 


3.74 


175 


10 


3.96 


2.28 




7 


2.45 


1.30 




14 


6.43 


3.96 


172.5 


10 


4.19 


2.42 




l-r 

( 


2.58 


1.38 




14 


6.76 


4.18 


170 


10 


4.42 


2.56 




7 


2.73 


1.47 



Table II: Inclusive single top quark production cross sections for top quark production (left) and 
antitop quark production (right) at E c . m . = 14, 10 and 7 TeV and three different m t . 

(a) (b) 

Top a. 4 Antitop 

— m T = 175GeV /,..-*• e" —1% = 175 GeV 

m T = 172.5 GeV / ■■■■ rrw = 172.5 GeV 

1% = 170 GeV 1% = 170 GeV 
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Figure 2: E c/rn . dependence of the single top (a) and antitop (b) quark production cross section for 
m T = 175, 172.5, 170 GeV. 

and E c . m = 900 GeV, the inclusive cross section for top (antitop) quark production is 0.018 
(0.005) pb. 

For m t = 175 GeV, Fig. |3] compares the E c , m , dependence of NLO and LO cross sections, 
for top and antitop quark, and for the charge asymmetry ratio (TZ). 1Z is defined as the 
ratio of top over antitop quark production cross sections. Different from the Tevatron, 
the single top quark production at the LHC (a pp collider) has a large charge asymmetry, 
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Figure 3: E c/m , dependence of single top (a) and antitop (b) quark production cross section at 
NLO and leading order (LO) for my = 175 GeV and ratio of cross sections for top/antitop quark 
production (c). 

leading to a difference in the numbers of top versus antitop quarks produced. Such an 
asymmetry is preserved in the charge of charged leptons from the top quark decay and will 
be measured at the LHC. It can be seen that for both top and antitop quark, with rising 
E cm , the cross section grows faster at NLO than at LO. This is mainly due to the higher 
momentum of the sea quarks and gluons in the colliding protons which increases the initial 
state corrections. The ratio of top/antitop quark production decreases with higher energy 
as the PDF difference between up and down quark looses significance. 

Without losing generality, in the remaining paper we consider E crai _ = 14 TeV only. 



B. Top quark mass dependence and theoretical uncertainties 

In order to predict the cross section for single top quark production as precisely as pos- 
sible, we need to understand how it depends on variations of the input parameters, such as 
m t dependence, scale dependences, and PDF uncertainties. 

The Tevatron has accomplished to reduce the uncertainty of the newest world average m t 
to 1.3 GeV jf)?]], so it is of interest to see how the cross section varies in this range. It can 
be seen in Fig. HJ that a variation of m t at 175 GeV of ± 5 GeV changes the cross section 
of single top quark production by about =p 10%. For the current world average of 173.1 ± 



1.3 GeV (57j, the predicted numerical result for the s-channel single top quark production 
cross section is 6.36 pb ± 0.19 pb, where the error of 3% is due to mass uncertainty only. 
For antitop quark production the calculation yields 3.90 pb ± 0.12 pb, also with an error of 
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Figure 4: nit dependence of single top (a) and antitop (b) quark production cross section and ratio 
of cross sections for top/antitop quark production (c), E cm _ = 14 TeV. 

3%. 

In Fig. H] (c) the charge asymmetry ratios are ~ 1.6 at LO and NLO in the region of 
168 GeV < m t < 182 GeV. Such a large charge asymmetry will be detectable at the LHC. 
We also note that the ratio curves at LO and NLO both rise with increasing m t , which can 
be understood as follows. If we set X\ ~ 22 = x where x± t 2 is the fraction of incoming proton 
energy carried by the parton, we obtain the average x as (x) ~ m t / 'y/s from s = X1X2S, 
where a/s = E c rn and y/l is the invariant mass of the incoming partons. The heavier 
the top quark, the larger (x). Since the down quark PDF peaks at lower x value than 
the up quark PDF, the antitop quark production rate decreases faster than the top quark 
production rate with increasing x (i.e. increasing m t ). Therefore, it yields the increasing 
ratio curves in Fig. |4](c). The 0(a s ) corrections, involving more production channels, only 
distort this picture slightly. For the remainder of this paper we use m t = 175 GeV. 

Besides m t dependence, single top quark production also suffers from scale dependence, a 
theoretical uncertainty originated from the unknown higher order corrections. There are two 
kinds of uncertainties: One is the renormalization scale fiR which is used for redefining the 
bare parameters in terms of the renormalized parameters, the other is the factorization scale 
fip which is introduced in order to absorb the collinear divergence into the PDFs. Although 
both fiR and fip are only introduced for technical reasons and our predictions for the cross 
section should not depend on them, we see such dependences, as we only work at first order 
in perturbation theory. In principle, /i# and //p are two independent theoretical parameters. 
For simplicity, we choose /i# = hf = Ho = m t . We vary /i by a factor 2 to estimate the 
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Figure 5: Scale dependence of single top (a) and antitop (b) quark production cross section and 
ratio of cross sections for top/antitop quark production (c). £ c . m . = 14 TeV and rrit = 175 GeV. 

size of higher order quantum corrections. Figure [5] (a, b) shows that the 0(a s ) corrections 
reduce the scale dependence of the single top and antitop quark production cross section. In 
Fig. |5] (c) we see that the charge asymmetry ratio falls at LO with rising fip, but increases 
at NLO. This difference comes from additional constituent processes available at NLO. 

Another theoretical uncertainty comes from the PDFs. The uncertainty of the 
CTEQ6.6M PDFs is given in Ref. 51] . At NLO the PDF uncertainty is found to be roughly 
about 3% for both top and antitop quark productions. Such a small uncertainty is due to rel- 
atively large x values (x ~ m t / 'y/s ~ 0.01) typical for s-channel single top quark production, 
where PDF uncertainties are generally small 51]. 



III. SINGLE TOP QUARK ACCEPTANCE STUDIES 

The W boson from the top quark can decay into jets or leptons. We only consider the 
leptonic decay mode as the hadronic decay mode is very difficult to observe experimentally, 
due to the large QCD background. Thus, the experimental signature of an s-channel single 
top quark event at NLO is the following: one charged lepton, missing transverse energy ($t) 
and two or three jets. As we are discussing single top quark production at parton level here, 
we only approximate the kinematic acceptance of the detector, and do not consider other 
detector effects such as ^-tagging efficiency or jet energy resolution. For the discussion of 
the effects of gluon radiation, jets have to be defined as infrared-safe observables. For this 



study we use the cone-jet algorithm [58j, as explained in Ref. [40| with cone-size R = 0.4 



59[. In this section we furthermore consider R = 1.0 for reference. The same i?-separation 
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2 & 3 jet 


1.42 


2.02 


0.57 


0.19 


-0.16 


0.90 


1.28 


0.37 


0.12 


-0.10 




3 jet 




0.81 


0.62 


0.15 


0.04 




0.50 


0.38 


0.09 


0.02 


(b) 


2 & 3 jet 


1.19 


1.69 


0.45 


0.17 


-0.12 


0.75 


1.08 


0.29 


0.11 


-0.07 




3 jet 




0.44 


0.38 


0.06 


0.01 




0.27 


0.23 


0.04 


0.004 


(c) 


2 & 3 jet 


0.68 


0.79 


0.21 


0.02 


-0.13 


0.44 


0.51 


0.14 


0.01 


-0.08 




3 jet 




0.24 


0.19 


0.05 


0.004 




0.15 


0.12 


0.03 


0.002 



Table III: LO and NLO contributions of cross sections for inclusive two-jet events as well as for 
exclusive three-jet events, for different cut scenarios, (a) 'loose' cuts with R cu t = 0.4, (b) 'loose' 
cuts with R cu t = 1-0, (c) 'tight' cuts with R cut = 0.4. E c . m . = 14 TeV and m t = 175 GeV. 

is also applied to the separation between leptons and jets. 

We consider two sets of kinematic cuts on the final state objects, a 'loose' and a 'tight' 



cut set [601 ]: 



p e T > 30 GeV , M < 2.5, 
$ T > 20 GeV , 

^>PT m , h\<v7 ax , (!) 

ARej > R C ut j ARjj > R cu t- 

For both sets of cuts, each event is required to have one lepton and at least 2 jets passing 
all selection criteria. The cut on the separation is chosen to be R cu t = 0.4 (and for reference 
1.0). The 'loose' set of cuts requires the jets to have a transverse momentum of at least 
pmm _ 3Q Q e y ; anc [ a pseudo-rapidity of at most rf^ ax = 5. At least one of the jets has 
to come from a b quark. The 'tight' set requires two b jets and all jets to have at least 
p ™* = 50 GeV and at most rj™ 1 * = 2.5. 

Table II III shows the cross sections for top and antitop quark production in the s-channel 
single top quark mode, split up into the different LO and NLO contributions after applying 
the two sets of cuts. For the 'loose' cut set the results are shown for two different values of 
Rcut and for the 'tight' cut set for R cut = 0.4. The results for the 'tight' cuts and R cu t = 1-0 



10 



follow the same trend as the other numbers in the table. For inclusive two-jet events and 
Rcut = 0.4, the top quark acceptance for the 'loose' cut set is around 33% both at LO and 
NLO level, and for the 'tight' cut set around 15% for the LO contribution and around 13% 
for the NLO contributions. For inclusive two-jet events and R cut = 0.4, the acceptance for 
antitop quarks is a bit higher than for top quarks: 34% at LO and NLO for the 'loose' cut 
set and 16% (14%) for the 'tight' cut set at LO (NLO). The low acceptance is mainly due 
to the pt cuts on the leptons and jets. The lepton is a product of the W boson decay and 
its pt distribution peaks under 30 GeV, cf. Fig. [7J Imposing a smaller p T cut on the lepton 
alone would not improve the overall acceptance much because the pt distribution of the b 
jet produced in association with the (anti)top quark also peaks around 30 GeV (cf. Fig. [H]). 
Hence, a low acceptance still follows from the p? cut on the b jet. A larger value for R cut 
reduces the cross section because more events fail the lepton-jet separation cut. For inclusive 
two-jet events after applying the 'loose' cuts, the reduction from R cut = 0.4 to R cut = 1.0 is 
about 17% for top and antitop quark production, both at LO and NLO level. We further 
note that the acceptance for antitop quark production is slightly higher than the acceptance 
for top quark production. This is due to the 77 cuts on the lepton as the lepton in antitop 
quark production remains in a more central rapidity region than the lepton in top quark 
production, cf. Fig. [JJ This follows directly from the fact that the top quark is boosted 
more than the antitop quark. The top quark, when produced, receives large contributions 
from the valence up quark in the large x regime, therefore it has large momentum along the 
moving direction of the incoming up quark. Such a boost effect is transferred to the top 
quark decay products, yielding the wide distribution of the lepton r\. 

Figure |6] (a, b) shows the NLO single (anti)top quark production cross section after 
applying the 'loose' cut set, for different rjj 10,20 as a function of the jet px threshold. It can 
be seen that an rf^ ax of 2.5 as applied in our 'tight' cut set does not reduce the acceptance 
much compared to the ^J 10 ^ = 5 of the 'loose' cuts, but the reduction to very central regions 
of the detector greatly decreases the acceptance. The fraction of three-jet events is relatively 
large even up to high jet px cuts, cf. Fig. |6] (c, d). Due to collinear enhancement from INIT 
corrections, the third jet has large pseudo-rapidity, so that the fraction of events with an 
additional jet greatly decreases when only considering central regions of the detector. In the 
following discussion of event distributions we will use the 'loose' cut set and R cut = 0.4. For 
the 'tight' cut set, the results are similar. 
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(a) (b) 
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Figure 6: Single top and antitop quark production cross section and fraction of three-jet events 
at NLO for varying jet pr cuts, after applying the 'loose' cut set. (a, b) Total cross section for 
inclusive two-jet events as a function of the jet pr cut for three different jet r\ cuts, (c, d) Fraction 
of exclusive three-jet events as a function of the jet pt cut for three different jet rj cuts. 

IV. SINGLE TOP QUARK EVENT DISTRIBUTIONS 

In this section we discuss some of the kinematic properties of final state objects in s- 
channel single top quark events. As described in Sec. IIII[ the final state consists of one 
charged lepton, fa, two b jets and possibly one additional jet due to NLO corrections. It is 
experimentally impossible to determine which of the b jets is produced in association with 
the top quark and which comes from the top quark decay. With the third jet from NLO 
corrections, the situation becomes even more complicated, as the additional jet has to be 
correctly identified to come either from the production or decay of the top quark. A third 
jet coming from the SDEC correction should be included in the top quark reconstruction, 
while a third jet from INIT or FINAL corrections should not be. To find the best prescrip- 
tion for classifying the gluon/light quark jet correctly, we first examine various kinematical 
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distributions of the final state particles. This includes a discussion of bfi n , the b (b) jet that 
is produced in association with the top (antitop) quark, and bd ec , which is the b (b) jet from 
the top (antitop) quark decay. We then describe different prescriptions for the top quark 
reconstruction. After showing the superiority of the best-jet algorithm, we are able to re- 
construct the top quark and discuss its kinematic properties. The kinematical distributions 
of the final state particles are examined after applying the 'loose' set of selection cuts as 
discussed in Eq. (JTJ). 

A. Final state object distributions 

1. Leptons and missing transverse energy 

Figure [7] shows a comparison between the kinematical distributions of lepton and fix 
for both top and antitop quark production at LO and with 0(a s ) corrections. As the 
distributions are for leptons and not quarks, the 0(a s ) corrections do not change the general 
form of the px distributions. The shapes of the distributions for the lepton p? look similar 
for top and antitop quarks. Furthermore, the $t distributions peak at about 32 GeV for 
both quark types. This is mainly because of a similar phase space in the transverse direction 
and similar spin correlations. 

It can also be seen that the $t distribution peaks at a higher energy than the lepton pt 
distribution. This is due to spin correlations and the left-handedness of the weak interaction 
and illustrated in Fig.[HJ The (anti) neutrino from the W boson decays preferentially following 
the moving direction of the (anti)top quark and is therefore boosted in comparison to the 
lepton, which decays in the opposite direction when seen in the rest frame of the W boson. 
The lepton 7] distributions are the same for positive and negative values of r\ as the pp initial 
state of the LHC is parity symmetric. The NLO 77 distribution for the positron from the top 
quark decay peaks at around rj = 1.4. This non-zero peak is due to LO and cannot be seen 
in the 0(a s ) corrections. It originates from the longitudinal boost that the intermediate 
W boson (Wi n t) receives from the PDFs. For example, a single top quark is preferentially 
produced through collision of an up valence quark and a down antiquark from the quark sea. 
The valence quarks carry a large momentum fraction of the incoming proton, while the sea 
quarks carry a small fraction. As a result, the (tb) system, which is equivalent to the (ud) 
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Figure 7: px (a, b) and rj (c, d) for the lepton and $t (e, f) after selection cuts, comparing LO 
to 0(a s ) corrections, (a, c, e) show top quark events, (b, d, f) antitop quark events. In (a, b) the 
distributions are also shown before the lepton pt cut. 

system at LO, is naturally boosted along the direction of the incoming valence quark. This 
longitudinal boost is less strong for antitop quark production as in the antitop quark case 
a down valence quark collides with an up antiquark, with the down valence quark carrying 
a smaller momentum fraction than the up valence quark. At NLO, the reacting parton in 
the initial state could be a sea (anti-) quark or gluon, which results in a smaller difference in 
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antitov moving direction 

Figure 8: Spin correlations for top (a) and antitop (b) quark decay in the (anti)top quark rest 
frame. The thin lines describe momenta while bold lines indicate spins. 

momenta between the colliding partons. Therefore, the lepton r\ distributions of the 0(a s ) 
corrections peak around zero. 

2. bfi n and b dec jets 

The comparison of the distributions of the bd ec and bf in jets is shown in Fig. [9] for 
single top quark production. It is experimentally not possible to distinguish b and b jets, 
but very instructive to consider their distributions individually. The px Born distribution 
of the bdec jet peaks at roughly 1/3 of m f , as it is a top quark decay product, while that of 
the bfi n jet peaks at lower p? and has a long tail to high p? values, because it is produced 
in association with the heavy top quark and has to balance the top quark p?. The peak 
positions of both bfi n and bd ec jets are shifted to a slightly lower value by the QCD corrections. 
This is mainly due to the FINAL (SDEC) contribution because the emitted gluon tends to 
be collinear to the bfi n (bd ec ) jet in the FINAL (SDEC) correction. Furthermore, the NLO 
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Figure 9: pt distributions for the bd ec jet from the top quark decay (a, c) and the bji n jet (b, d) 
after selection cuts [(b) is also shown before the bfi n jet pt cut], (a, b) show NLO, Born and sum 
of 0{a s ) contributions, (c, d) individual 0(a s ) contributions. 



distribution is broadened due to the INIT contribution which adds additional px to the event 
by emitting a third jet in the initial state. Note that among the three 0(a s ) corrections 
the INIT contribution dominates as it receives both soft and collinear enhancements at 
NLO. The soft gluon contributions in the FINAL and SDEC contributions are suppressed 
in comparison to the INIT correction, due to the large top quark mass. The single antitop 
quark production exhibits very similar distributions and is not shown here. 

Figure [10] shows the r] distributions of the bf in and bd ec jet for both single top and anti- 
top quark production. For top quark production the Born distribution of the bdec jet t] is 
much broader than the Born distribution of the bfi n jet t]. On the other hand, for antitop 
production the bfi n jet exhibits a wider 77 distribution than the bd ec jet. This is due to spin 
correlations and the boost of Wi n t along the beam line. Figure [UJ shows single top and 
antitop quark production in the cm. frame of W int . Due to spin conservation, for top quark 
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Figure 11: Pictorial illustration of spin correlations and boost effects in top (a) and antitop (b) 
quark production in the cm. frame of Wi n t- The thin lines denote particle momenta while the 
bold lines label spin directions. 
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Figure 12: Invariant mass of the (bdec jet, 6/« n jet) system after selection cuts, comparing Born 
level to 0(a s ) corrections. 

production, in the cm. frame of Wi n t, the top quark mainly follows the direction of motion 
of the incoming up quark, while the bfi n jet \b quark in Fig. [11] (a)] follows the direction 
of motion of the incoming antiquark. Furthermore, in the top quark production process, b 
is right-handed and b is left-handed due to the left-handed charged current interaction in 
the SM. Since the valence quarks carry a larger fraction of the proton momentum, Wi n t is 
boosted along the direction of motion of the up quark. The top quark and its decay product, 
the bdec jet, go in the same direction, and therefore have a larger longitudinal momentum and 
higher r\ than the bf in jet, which, in the cm. frame of Wi nt , follows the opposite direction 
and therefore receives a smaller longitudinal boost. For the antitop quark production it is 
the bfi n jet [b quark in Fig. [TT] (b)] and not the antitop quark that, in the cm. frame of 
Wi n t, follows the direction of motion of the incoming quark. Therefore the rj distribution of 
the bfi n jet is wider than the same distribution of the bdec jet from the antitop quark decay. 
As mentioned in Sec. IIV A 11 the boost is smaller for antitop quark production, so the effect 
on the t] width is less strong. 

The 0(a s ) corrections shift the bfi n jet produced in association with the top quark to more 
central rapidities. This effect is not very strong and mainly due to initial state corrections, 
which add to the pt of the event while reducing the boost along the beam line direction. 
The shape of the r\ distribution of the bd ec jet remains practically unchanged by the 0(a s ) 
corrections. 

Single top quark production is an irreducible background to the Higgs boson search, e.g. 
W ± H associated production with the subsequent Higgs boson decay H — > bb. Even though 
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W ± H production is not the largest production channel at the LHC, one has to combine 
its contribution with other channels in order to reach 5<r statistics significance. Thus it is 
crucial to have a good understanding of the SM backgrounds. Since a light Higgs boson 
predominantly decays into a bb pair, we examine the impact of 0(a s ) corrections on the 
invariant mass distribution of the bb pair in single top quark production. Figure [T21 plots the 
invariant mass distribution of the (bdec jet, b^ n jet) system. It can be seen that the FINAL 
correction shifts the peak of the invariant mass to slightly lower values. This is the case if a 
third jet is produced in addition to the bd ec jet and the bf in jet. The INIT corrections tend 
to have peaks at higher invariant mass, due to the additional pt they provide. The drop-off 
at higher mass values is faster at NLO level. The invariant mass distribution of the two 
b jets for antitop quark events looks very similar in spite of the differences in initial state 
PDFs. 



3. Invariant mass M(bd ec jet, lepton) 

The top quark decays into a bd ec jet and a W boson, which itself can decay leptonically. 
The invariant mass distribution of the (bd ec jet, lepton) system is characteristic of the d ecay 



of the SM top quark and sensitive to the W-t-b coupling (or the W boson helicity) [6l|, |62| , 



m b jet e ~ 2 ( m t ~ m w) ( X - cos ^) 

where d\ is the polar angle of the charged lepton in the rest frame of the W boson which 
is defined in the rest frame of the top quark. Figure [13] shows the drop-off behavior of the 
invariant mass distribution at x/mf — m 155 GeV, due to the kinematics of the event. 
It also shows that the SDEC correction shifts the invariant mass peak to lower values, mainly 
because it weakens the spin correlations. 



4- Ht distributions 

In order to distinguish the s-channel single top quark mode from dominant backgrounds, 
it is important to know how the total transverse energy Ht changes with NLO corrections. 
In Fig. [14] we look at the effect that 0(a s ) corrections have on the total transverse energy 
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(a) (b) 




H T [GeV] H T [GeV] 

Figure 14: Comparison of the total transverse energy of single top quark events after selection cuts 
between LO and NLO (a) and between the different 0(a s ) corrections (b). 



Ht of a single top quark event. Ht is denned as 

H T =p l e pton +# T + YsPT et - (2) 

jets 

Clearly, the SDEC and FINAL 0(a s ) contributions shift the total transverse energy down, 
while the INIT contribution shifts it up. This again is due to the additional transverse energy 
that a third jet in the initial state adds to the event. The shift to higher transverse energies 
is significant for single top quark measurements at the LHC, as tt events are generally a 
larger background for higher Ht values. 
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B. Event reconstruction 



One of the main reasons for studying single top quark events is to find out more about the 
properties of the top quark and its couplings. Furthermore, as single top quark production 
is a weak interaction process, we expect a number of correlations between the particles. 
It is therefore of interest to reconstruct the complete event, including final state jets and 
intermediate particles. 



As explained in Refs. 40|, |63], it is possible to reconstruct the W boson from the observed 



charged lepton and fix, where the unknown longitudinal momentum of the neutrino p v z is 
substituted by the restriction that the invariant mass of the (charged lepton, neutrino) 
system has to equal the mass of the W boson. Of the two possible solutions we pick the 
one with the smaller \p v z \. For the top quark reconstruction, we must then combine the 
reconstructed W boson with the bd ec jet from the top quark decay. This means we have to 
identify the correct jet as the bd ec jet. There are several different methods to select this jet. 
As both the bdec jet and the bfi n jet have high pt and are possibly b-tagged, neither the 
method of choosing the leading jet (jet with the highest px) nor the second jet (jet with the 
second highest p?) nor the b-tagged jet is very reliable in identifying the right jet. 

A more effective algorithm is the so-called best-jet algorithm as explained in Ref. 64] 



and more specifically for our case in Ref. 



40| | . Here, the Wj or the Wjj (only possible in 



cases where there is a third jet) combination, that gives an invariant mass closest to the 
input m t , is chosen as the reconstructed top quark and the jet (or two-jet system), that is 
thus identified, is called best-jet. Of the two possible other jets, the one with larger p T is 
the so-called non-best-jet. 

Figure [T5l shows the efficiencies of the best-jet, the leading jet (jet 1) and the second jet 
(jet 2) algorithms. The leading jet has an overall efficiency of 39%. The second jet mainly 
corresponds to the bd ec jet for very high p^ and has an overall efficiency of about 49%, 
while the best-jet algorithm has a high efficiency for all momenta and identifies the bd ec jet 
correctly in 93% of all events. Its effectiveness is mostly limited by the efficiency of the W 
boson identification; with a falsely reconstructed W boson, the identification of the bd ec jet 
by the best-jet algorithm becomes a random pick. 

The invariant mass of the Wj combination is shown in Fig. [161 comparing the reconstruc- 
tion using the best-jet, the leading jet and the second jet algorithm to the reconstruction 
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b decJ et P T [ GeV ] 



Figure 15: px distribution of the bdec jet from the top quark decay, for all events that pass the 
'loose' set of cuts (solid), only for those events in which the bdec jet is the best-jet (dashed), only 
for events in which the bdec jet is also the leading jet (dotted) and only for events in which the bdec 
jet is the second jet (dash dotted). 
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Figure 16: Invariant mass of the Wj combination after selection cuts, where the jet j is the real 
bdec jet (solid), the best-jet (dashed), the leading jet (dotted) and the second jet (dash dotted). 

using the real bd ec jet. For the identification of the real bd ec jet, parton level information is 
used and a possible third jet from the top quark decay is included. For all four curves, the 
W boson is reconstructed from the final state lepton and $t, so that all differences in shape 
and height are due to the method of identifying the bdec jet. As expected, the distribution 
using the best-jet fits the distribution with the true bd ec jet information much better than 
the distribution using the leading jet or the second jet. It is important to notice though, 
that this is due to two competing effects: (i) The best-jet does not identify the correct bd ec 
jet in 7% of all events, which reduces the height of the distribution in Fig. [15] (ii) Due to the 
requirement to be as close as possible to the input m t , mis-reconstructed events are shifted 
closer to 175 GeV, increasing the height of the peak. The second effect is larger than the 
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Figure 17: Distributions of px (a, b) and rapidity (c, d) of the top quark (a, c) and of the antitop 
quark (b, d), both reconstructed with the best-jet algorithm and after applying the 'loose' set of 
selection cuts. 



first, so that the peak reconstructed with the best-jet algorithm is actually higher than the 
peak reconstructed using the real bd ec jet. 

We reconstruct the top quark using the best-jet algorithm and are now able to study 
some of its kinematic properties. In Fig. [T71 the px and rapidity distributions are shown for 
both top and antitop quark, reconstructed with the best-jet algorithm and after applying the 
'loose' set of selection cuts. The shapes of the pt distributions look similar, but the rapidity 
distribution is wider for the top quark than for the antitop quark. This is a reflection of the 
boosted kinematics of Wi nt and the (anti)top quark itself (cf. Fig. ITT]) . The NLO corrections 
do not significantly change the the shape of the distributions. 

We reconstruct the virtual W boson Wi nt in single top quark production by combining 
the reconstructed top quark with the non-best-jet. This method is exact if best-jet and 
non-best-jet are identified correctly and the event does not contain a third jet from 0(a s ) 
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Figure 18: p z (a, b), rapidity (c, d) and invariant mass (e, f) distributions of the intermediate W bo- 
son after selection cuts, (a, c, e) show top quark production, (b, d, f ) antitop quark production. The 
rapidity distribution only contains events with px ^ 0. 



FINAL corrections. As discussed in Sec. IIV A 11 at LO Wi n t receives a longitudinal boost, 
which is stronger for top than for antitop quark production. Figure [18] shows that the LO 
p z and rapidity distributions are wider for Wi nt in top than in antitop quark production. 
Furthermore, the 0(a s ) correction to the rapidity distribution of Wi nt is narrower than the 
LO contribution. This is because at NLO the incoming partons could also be sea quarks or 



24 



(a) 



(b) 



> 10 F Top 

C5 -NLO 

5 8 Born 

JO I 0(a s ) sum 



Q. 

O 




-200 -100 



o 

JO 

< 
e 



100 
[GeV] 



12 
10 

8 
6 
4 
2 





Top 




-NLO 




Born 




0(ci s ) sum 





i i i""; i-1 ... 



12 3 4 

Kiet J - T|(jet_ 



Figure 19: Difference in px (a) and r\ (b) between the best-jet and the non-best-jet, after applying 
the 'loose' set of kinematic cuts. 



gluons which results in a less boosted Wi nt system. 

In order to search for new physics in the form of a W boson, it is important to know 
the invariant mass distribution of the intermediate SM W boson. W boson searches in 
the single top quark final state have been performed at the Tevatron and have set lower 
W boson mass limits in the range 750 GeV to 800 GeV 18l-l20l|. The LHC with its higher 
E c .m. will allow to probe up to even larger boson masses. The invariant mass distribution 
of Wi n t is shown in Fig. [18] for SM single top and single antitop quark production. The 
0(a s ) corrections, more specifically the FINAL contribution, shift the invariant mass peak 
to lower energies, as in our reconstruction of Wi nt we do not include a possible third jet in 
the final state. 



C. Kinematical and spin correlations 

1. Correlations between the b jets 

Having reconstructed and identified the W boson, the b^ec jet and the bf in jet, it is now 
possible to study correlations between objects in single top quark events. Figure [19] shows 
that the bd ec jet (best-jet) and the bfi n jet (non-best-jet) in single top quark events are highly 
correlated. The correlations are modified by the 0(a s ) corrections. The FINAL corrections 
shift the pt difference between best-jet and non-best-jet to more positive values, as their 
tendency to be collinear to the boost direction lowers the pt of the bfi n jet. For the same 
reason, the r] correlation is shifted to slightly more central values. 
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top quark production top quark decay 




Figure 20: S-channel single top quark production from the perspective of a top quark at rest. Thin 
arrows denote directions of momentum while bold arrows indicate spin directions. 

2. Top quark polarization 

The SM predicts that single top quarks are highly polarized. Verifying this prediction 
can be used as a test of the SM electroweak symmetry breaking mechanism and a check 



for new physics 



23 



24 



651 ] . It is in principle possible to measure the polarization of the 



single top quark, by making use of the fact that the charged lepton from the top quark 



decay is maximally correlated with the top quark spin |66l . |67| , as illustrated in Fig. [HJ In 
the following we therefore plot the angle between the charged lepton and a reference axis 
in the rest frame of the top quark. Two different choices for this reference frame have been 
used in the past: the helicity basis and the so-called optimal basis jsl, In the helicity 

basis, the top quark spin is measured along the top quark direction of motion in the center 
of mass (cm.) frame of the system. In the optimal basis, the top (antitop) quark spin is 
measured along the direction of the incoming antiquark (quark) in the cm. frame of the top 
(antitop) quark. This reference frame is called the optimal basis, as the top (antitop) quark 
produced in s-channel single top quark processes is almost 100% polarized along (against) 
the direction of the incoming antiquark (quark). Figure [201 shows s-channel single top quark 
production from the perspective of a top quark at rest. It demonstrates how the optimal 
basis makes use of the momentum - spin correlations that are due to spin conservation and 
left-handedness of the weak interaction. 

While at the Tevatron the antiquark comes predominantly from the antiproton, the LHC 
is a pp collider and there is no preferred direction for the antiquark. We can however make 
use of the fact that the incoming quark is more likely to be a valence quark, and has a larger 
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longitudinal momentum than the antiquark which comes always from the quark sea of the 
proton. In most cases the direction of the longitudinal boost of Wi n t therefore indicates the 
direction of the incoming quark. (This becomes a better approximation when the magnitude 
of the longitudinal momentum \p z \ of Wi nt is large.) We choose the reference axis according 
to the sign of p z {Wi nt ). After discussing the helicity basis, we show how it is possible to 
enhance the spin correlations for the optimal basis further by requiring a minimum p z (Wi n t). 

For the helicity basis, the cm. frame of the system has to be reconstructed in order 
to define the top quark momentum. As discussed in Ref. [41 , this is more complicated at 
NLO than at Born level, because of additional jets. Ref. [41] investigates two options for 
reconstructing the cm. frame: the tbfi n (j)-frame, which is the cm. frame of the incoming 
partons and the rest frame of all the final state objects (reconstructed top quark and all other 
jets), and the £&/; n -frame, which is the cm. frame of the top quark and the non-best-jet. 
The tbfi n -fr&me differs from the tbf in (j)-iiame only in exclusive three-jet events. As shown 
in Table fVl and discussed below, the degree of polarization is larger in the t&/; n -frame, so we 
only show distributions for the top quark polarization in this frame. 

In the helicity basis, we examine the polarization of the top quark by studying the angular 
distribution cos 9hei of the lepton relative to the moving direction of the top quark, both in 
the cm. frame of the system, 



cosO hel = — — -, (3) 

\pt\W 



Pt _ 

\pi\ 

where pt is the top quark three-momentum defined in the tb^ n - or the t6/j n (j')-frame, and 
pi is the charged lepton three-momentum, after boosting it first into the cm. frame of the 
system and then into the top quark rest frame. For a left-handed top quark, the angular 
correlation of the lepton £ + is given by (1 — cos6 l / ie ;)/2, and for a right-handed top quark, it 
is (1 + cos6^ e ;)/2. For a right-handed antitop quark, the angular correlation of the lepton 
l~ is given by (1 — cos 9hei)/2, and for a left-handed antitop quark, it is (1 + cos Ohei)/^- 

Figure [21] shows that this linear relationship for cos &hei indeed describes s-channel single 
top quark events well at parton level. It also shows that the top quark is not completely 
polarized in the helicity basis, and that this polarization is weakened further at NLO. After 
event reconstruction, the drop-off close to cos Ohei = -1 is due to jet-lepton separation cuts. 

The corresponding distributions for the optimal basis can be seen in Fig. [221 Here, the 
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Figure 21: Top quark polarization in the helicity basis at parton level (a, b) and after event 
reconstruction and applying the 'loose' set of selection cuts (c, d), for top quark production (a, c) 
and antitop quark production (b, d). The tb /-^-frame is chosen as the cm. frame. 



relevant angular correlation is 



cos 6, 



opt 



Ppl - Pi 



PI 



(4) 



where p Pl is the three-momentum in the top quark rest frame of the proton which travels (in 
the case of top quark production) antiparallel to the longitudinal boost that Wi nt receives, 
e.g. in the — pl(Wi nt ) direction. In the case of antitop quark production the chosen proton 
travels parallel to the longitudinal boost of W int , e.g. in the p" z (W int ) direction. p| is 
the three-momentum of the lepton in the top quark rest frame. For a(n) (anti)top quark 
polarized along (against) the moving direction of the chosen proton, the angular distribution 
of the lepton £ + (£~) is (1 + cos# OJrf )/2, while for a(n) (anti)top quark polarized against 
(along) the moving direction of the chosen proton, it is (1 — cos6 opt )/2. Figure I2"2l shows 
that there is indeed a linear relationship for cos 8 opt at parton level, but that the top quark 
is not completely polarized in this basis. After reconstruction and applying the 'loose' set 
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(a) (b) 




Figure 22: Top quark polarization in the optimal basis at parton level (a, b) and after event 
reconstruction and applying the 'loose' set of selection cuts (c, d), for top quark production (a, c) 
and antitop quark production (b, d). 

of cuts, there is a cutoff at large cos8 opt , due to the lepton 77 cut. 

It is possible to enhance the performance of the optimal basis, by using only those events 
in which p z (Wi n t) is larger than a certain threshold. For this study, thresholds of 500, 1000 
and 2500 GeV have been tested. By imposing a cut on p z (Wi n t), we are able to determine 
the direction of the incoming antiquark correctly more often, as the momentum difference 
between the incoming quarks increases. The efficiencies for the three different thresholds 
500, 1000 and 2500 GeV are 50%, 25% and 2% respectively (cf. Fig.UBJ. Figure [23] shows 
that, at parton level, the (anti)top quark is indeed highly polarized for a p z (Wi nt ) cut of 1000 
GeV. Reconstruction and selection cuts change the distribution dramatically and differently 
for top and antitop quarks. The cutoff close to cos 6 op t = 1 is again due to the lepton rj 
cut. Due to spin correlations, the lepton px cut carves out a large number of events with 
cos# opi < for the top quark, while for the antitop quark it removes mainly events with 
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Figure 23: Top (a, c) and antitop (b, d) quark polarization in the enhanced optimal basis, where 
Wi n t has a longitudinal momentum of at least 1000 GeV. (a, b) show the distributions before 
selection cuts using all parton information, while (c, d) are the distributions for reconstructed 
events after applying the 'loose' set of cuts. 

cos 6 ' pt > 0. 
As in Ref. 



4=0l ] . we define the degree of polarization T> of the top quark as the ratio 

AL — iV+ 



V 



N^ + N+ 

where iV_ (N + ) is the number of left-hand (right-hand) polarized top quarks in the helicity 
basis. Similarly, in the optimal basis, N_ (N + ) is the number of top quarks with polarization 
against (along) the direction of the chosen proton three momentum in the top quark rest 
frame p Pl . For the antitop quark the relationships are the same. 



The spin fractions T± and the asymmetry A of the distribution are defined in Ref. [40] . 
Without imposing any kinematic cuts, D = 2A, which can indeed be seen in Table IIVL 
Furthermore, the ratio of top quarks with spin along the basis direction will be r T = 0.5 — A 
when no cuts are applied. However, when cuts are imposed, the two relationships break 
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Top 


Ant it op 


V 

LO NLO 


T 

LO NLO 


A 

LO NLO 


V 

LO NLO 


T 

LO NLO 


A 

LO NLO 


Helicity (i&/j n -frame) 
Optimal (p z (W int ) > 0) 
Optimal ( Pz (W int ) > 500) 
Optimal (p z (W int ) > 1000) 
Optimal ( Pz (W int ) > 2500) 


0.69 0.68 
-0.61 -0.63 
-0.79 -0.82 
-0.88 -0.89 
-0.98 -0.97 


0.84 0.84 
0.81 0.81 
0.90 0.91 
0.94 0.95 
0.99 0.98 


0.34 0.34 
-0.31 -0.32 
-0.40 -0.41 
-0.44 -0.45 
-0.49 -0.48 


-0.68 -0.67 
0.53 0.52 
0.74 0.73 
0.83 0.81 
0.93 0.90 


0.84 0.83 
0.77 0.76 
0.87 0.86 
0.91 0.91 
0.96 0.95 


-0.34 -0.33 
0.27 0.26 
0.37 0.36 
0.41 0.41 
0.46 0.45 



Table IV: Parton level degree of polarization T>, polarization fraction J 7 , and asymmetry A for s- 
channel single top quark events before any cuts. Results are shown for both top (left) and antitop 
(right) quark measured in the helicity basis (t6jj n -frame) and in the optimal basis with different 
Pz(Wi n t) thresholds. In this table, T corresponds to T— in the helicity basis for left-handed top 
quarks and to T+ in the optimal basis for top quarks with polarization along the chosen proton 
three-momentum. For right-handed antitop quarks, T corresponds to JF+ in the helicity basis 
and to T- in the optimal basis for antitop quarks with polarization against the chosen proton 
three- momentum. 

down. 

Table [IV] shows T>, J 7 and A at parton level before any cuts. It can be seen that in 
the helicity basis the polarization is slightly larger for top than for antitop quarks, as top 
quarks receive larger longitudinal boosts from Wi nt and therefore have higher energies (cf. 
Sec. IIV Aip . The comparison between LO and NLO shows a slight decrease in \T>\ in the 
helicity basis for both top and antitop quark. In the optimal basis, the larger longitudinal 
boost also increases the energy of the top quark, but more importantly makes it more likely 
that the chosen proton is the correct pick for the reference axis. As shown in Table IIV] 
both quark types are more polarized in the helicity basis than in the optimal basis, if no p z 
cuts on Wi n t are applied. With higher p z cuts on W int , \T>\ calculated in the optimal basis 
increases though, as higher cuts improve the chance of guessing the correct direction of the 
initial antiquark. 

The results after reconstructing the jets {jP T > 30 GeV, \rjj\ < 5, ARej = ARjj > 0.4) 
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T O NT D 


T O 


NT O 


Helicity (tbf in (j), incl. 2-jet) 


0.74 


0.67 


U.o/ U.o4 


0.37 


0.34 


-0.74 


-0.54 


U.o f U. II 


-0.37 


-0.27 


Helicity (tbfi n , incl. 2-jet) 


0.74 


0.72 


0.87 0.86 


0.37 


0.36 


-0.74 


-0.71 


0.87 0.86 


-0.37 


-0.35 


Helicity (tb fin (j), excl. 3-jet) 




0.76 


0.88 




0.38 




-0.65 


0.83 




-0.33 


Helicity (tbf in , excl. 3-jet) 




0.81 


0.91 




0.41 




-0.80 


0.90 




-0.40 


Optimal (incl. 2-jet) 


-0.66 


-0.67 


0.83 0.83 


-0.33 


-0.33 


0.59 


0.56 


0.79 0.78 


0.29 


0.28 


Optimal (excl. 3-jet) 




-0.69 


0.85 




-0.35 




0.58 


0.79 




0.29 



Table V: Degree of polarization T>, polarization fraction J 7 , and asymmetry A for inclusive two-jet 
and exclusive three-jet s-channel single top quark events after jet clustering (p 3 T > 30 GeV, \r]j\ < 5, 
ARgj = ARjj > 0.4). Results are shown for both top (left) and antitop (right) quark measured in 
the helicity basis comparing the two different cm. frames, and in the optimal basis. In this table, 
T corresponds to T— in the helicity basis for left-handed top quarks and to JF+ in the optimal 
basis for top quarks with polarization along the chosen proton three-momentum. For right-handed 
antitop quarks, J- corresponds to in the helicity basis and to J-- in the optimal basis for antitop 
quarks with polarization against the chosen proton three-momentum. 

are listed in Table |V] for inclusive two-jet events as well as exclusive three-jet events. In 
comparison to the results before any cuts, \D\ is generally larger here, as the px cuts on the 
jets increase the energy of the (anti)top quark. Similarly, in both bases the polarization at 
NLO is higher for exclusive three-jet events than for inclusive two-jet events, because of the 
additional pr that a third jet in the initial state adds to the event, which again increases the 
energy of the (anti)top quark. Furthermore, the table shows that for the helicity basis the 
tbfi n -frame leads to higher polarizations than the tbfi n (j)-fiame at NLO. This is because in 
most of the three-jet events the third jet comes from INIT corrections and should not be 
included in the cm. system of the final state objects produced from Wi nt . 

Figure [2H shows the degree of polarization T> of top and antitop quark in the optimal 
basis, for different values of p z {Wi n t) at parton level for top and antitop quark production. 
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(a) (b) 




W int p 2 [GeV] W int p z [GeV] 

Figure 24: Degree of polarization T> of top quark (a) and antitop quark (b) in the optimal basis, 
for different values of p z (Wi n t), at parton level before cuts. 

D. Distributions for three-jet events 

Single top quark events at the LHC contain a large fraction of three-jet events, which 
can be seen in Fig. [6j It is therefore of interest to discuss the kinematic properties of the 
third jet. After selection cuts, this third jet corresponds to 0(a s ) corrections in about 80% 
of the three-jet events. This means, for most events, the emitted gluons/light quarks have a 
lower pt than the bd ec and bfi n jets, which is mainly due to the large amount of initial state 
radiation that tends to be collinear to the beamline. 

The emission of additional gluons/light quarks can be divided into production-stage emis- 
sion and decay-stage emission. Production-state emission includes INIT and FINAL cor- 
rections and occurs before the top quark goes on-shell, while decay-stage emission consists 
of the SDEC contribution and occurs after the top quark goes on-shell. This classification 
of three-jet events into production-stage or decay-stage is useful, but blurred by the finite 
width of the top quark and, in experiments, by jet-energy resolution and ambiguities in jet 
assignment. 

1. Kinematic distributions of the extra jet 

The pt and rj distributions of the third jet for different NLO corrections are presented 
in Fig. [251 The dominance of the initial state corrections mentioned above is due to the 
collinear enhancement of the incoming partons and can be seen in both distributions. As 
it is determined by the collider energy, the INIT contribution extends to far higher pt and 
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Figure 25: (a) and r/ distribution (b) of the third jet after applying the 'loose' set of cuts for 
the different 0(a s ) corrections. 



(a) 



(b) 



Top 

-0(a s ) sum 
INIT 
FINAL 
SDEC 




Top 

-0(a s ) sum 
INIT 
FINAL 
SDEC 



AR <i e4t 3'i e, bes ti J 




2 4 

AR(jet Jet 

u 3 ' non-best jet 



Figure 26: Difference in the pt between the third jet and the best-jet, if the best-jet is not the 
third jet (a) and the third jet and the non-best-jet (b), after applying the 'loose' set of kinematic 
cuts. 



larger rj than the other contributions. 

Figure [261 shows that it is in principle possible to identify the third jet as coming from the 
SDEC or the FINAL corrections, by looking at the difference in the px between the third jet 
and the best-jet or the third jet and the non-best-jet respectively. The pt difference between 
the third jet from the decay and the best-jet has a peak close to zero for the decay contribu- 
tion, while the difference between the px of the third jet from the FINAL correction and the 
non-best-jet tends to be smaller than for the other contributions. Unfortunately those peaks 
are not very high, and the initial state correction is again dominant for slightly larger pt 
differences, so that an optimal cone size has to carefully balance out the competing effects of 
desired jets falling outside of the cone and unwanted initial state radiation being included. 
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Figure 27: Angular correlation between the third jet and the best-jet after selection cuts. All NLO 
corrections except for SDEC (solid), denoted as 0(a s )nd, are compared to the SDEC contribution 
only (dotted). Here, we only allow events in which exactly two jets are used to form the best-jet. 

This is further complicated in experiments by hadronization and detector resolution effects. 



2. Angular correlation between the extra jet and the best-jet 

The best-jet algorithm can be used to distinguish production-stage from decay-stage 
parton (gluon or quark) emission. In Fig. [23 we show the angular correlation between the 
third jet and the best-jet after event reconstruction. Here, we include only events in which 
exactly two jets are used to form the best-jet. There is a clear separation between the SDEC 
contribution only and the rest of the NLO corrections, even after event reconstruction with 
kinematic cuts imposed. 



V. CONCLUSIONS 



We have presented a study of s-channel single top and antitop quark production at the 
LHC based on the full NLO calculations. We have studied the NLO QCD corrections to the 
production and decay of the single top quark and shown their effect on the inclusive cross 
section as well as on kinematical distributions, including top quark polarization measure- 
ments. For this, we have divided the higher order corrections into three Gauge invariant 
sets. The inclusion of NLO corrections allows more precise predictions of the properties of 
single top quark events, which is mandatory for using single top quark events to test the 
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SM and search for new physics. 

The NLO corrections increase the s-channel inclusive cross section significantly. Due to 
the PDF of the proton, the cross section is larger for single top quark production than for 
single antitop quark production and depends on the value of m t . Simple kinematic cuts, as 
they are used by the ATLAS and the CMS experiments, reduce the acceptance considerably, 
while the percentage of three-jet events remains high, especially due to collinear enhancement 
of the initial state corrections. 

The LO kinematical distributions of final and intermediate state objects in single top 
quark events are dominated by spin correlations and the momentum difference between the 
two incoming partons. These effects are generally smaller for antitop quark production than 
for top quark production and reduced by 0(a s ) corrections in both cases. This is because 
the NLO corrections weaken spin correlations and, in some cases, lower the momentum 
difference between the incoming partons. 

The total transverse energy Ht, the invariant mass of the {bdec jet, bf in jet) system, and 
the invariant mass of the (bd ec jet, lepton) system are examples of distributions that are 
characteristic of single top quark events. They change significantly if NLO corrections are 
included. The invariant mass distribution of the SM Wi nt could be a useful discriminator 
for W boson searches at the LHC. 

For the reconstruction of the top quark it is important to identify the correct jet as the 
bdec jet, that is the jet which is produced when the top quark decays. The most efficient 
method for this identification is found to be the best-jet algorithm, which picks the Wj or 
Wjj combination that gives an invariant mass closest to the input m t . 

The spin correlations and the fact that the top quark decays before it can hadronize makes 
it possible to measure its polarization. We identify appropriate frames for top quark spin 
correlation measurements and find that an additional cut on p z (Wi nt ) significantly increases 
the measured spin correlations in the optimal basis. 



36 



Acknowledgments 

S. H. and R. S. are supported in part by the U.S. National Science Foundation under 
Grant No. PHY-0757741. Q. H. C. is supported in part by the Argonne National Laboratory 
and University of Chicago Joint Theory Institute (JTI) Grant No. 03921-07-137, and by the 
U.S. Department of Energy under Grants No. DE-AC02-06CH11357 and No. DE-FG02- 
90ER40560. C. P. Y. acknowledges the support of the U.S. National Science Foundation 
under Grant No. PHY-0555545 and PHY-0855561. C. P. Y. would also like to thank the 
hospitality of National Center for Theoretical Sciences in Taiwan and Center for High Energy 
Physics, Peking University, in China, where part of this work was done. 



DO, V. M. Abazov et al, Phys. Rev. Lett. 98, 181802 (2007), arXiv:hep-ex/0612052 



DO, V. M. Abazov et al, Phys. Rev. D78, 012005 (2008). la7XTv:0803.0739l 

CDF, T. Aaltonen et al, Phys. Rev. Lett. 101, 252001 f2008b la7XTv:0809.25811 

CDF, T. Aaltonen et al, Phys. Rev. Lett. 103, 092002 f2009b larXTv:0903.0885l 

DO, V. M. Abazov et al, Phys. Rev. Lett. 103, 092001 f2009b larXTv:0903.0850l 

K.-i. Hikasa, K. Whisnant, J. M. Yang, and B.-L. Young, Phys. Rev. D58, 114003 (1998), 



arXiv:hep-ph/9806401 



T. Han, M. Hosch, K. Whisnant, B.-L. Young, and X. Zhang, Phys. Rev. D58, 073008 (1998), 



arXiv:hep-ph/9806486 



T. Tait and C.-P. Yuan, Phys. Rev. D63, 014018 (2001), [arXiv:hep-ph/0007298 

J. J. Liu, C. S. Li, L. L. Yang, and L. G. Jin, Nucl. Phys. B705, 3 (2005), 



arXiv:hep-ph/0404099 



O. Cakir and S. A. Cetin, J. Phys. G31, Nl (2005). 

J. Guasch, W. Hollik, S. Penaranda, and J. Sola, Nucl. Phys. Proc. Suppl. 157, 152 (2006), 



arXiv:hep-ph/0601218 



D. O. Carlson, E. Malkawi, and C.-P. Yuan, Phys. Lett. B337, 145 (1994), 



arXiv:hep-ph/9405277 



T. Tait and C.-P. Yuan, Phys. Rev. D55, 7300 (1997), [ailav:hep-ph/96lT244 



D. Espriu and J. Manzano, Phys. Rev. D65, 073005 (2002), arXiv:hep-ph/0107112 



37 



Q.-H. Cao, C. S. Li, and C. P. Yuan, Phys. Lett. B668, 24 (2008), [aTXiv:hep-ph/0612243 



Q.-H. Cao, J. Wudka, and C. P. Yuan, Phys. Lett. B658, 50 (2007), arXiv:0704. 2809, 
J. A. Aguilar-Saavedra, Nucl. Phys. B804, 160 (2008), larXiv:0803.381~0l 



DO, V. M. Abazov et al, Phys. Lett. B641, 423 (2006), |arXiv:hep-ex/0607102 



DO, V. M. Abazov et al, Phys. Rev. Lett. 100, 211803 (2008). la7XTv:0803.32561 



CDF, T. Aaltonen et al, (2009), arXiv:0902.3276 



DO, V. M. Abazov et al, Phys. Rev. Lett. 99, 191802 (2007), arXiv:hep-ex/0702005 



CDF, T. Aaltonen et al, Phys. Rev. Lett. 102, 151801 (2009). la7XTv:0812.34001 
DO, V. M. Abazov et al, Phys. Rev. Lett. 101, 221801 (2008). larXTv:0807. 16921 
DO, V. M. Abazov et al, Phys. Rev. Lett. 102, 092002 (2009). larXTv:0901.0151l 
DO, V. M. Abazov et al, Phys. Rev. Lett. 102, 191802 (2009). larXTv:0807.08591 



S. Moretti, Phys. Rev. D56, 7427 (1997), [aiXiv:hep-ph/970 5388 
K. Odagiri, (1999), [aTxTv:hep-ph/990l"432 



M. C. Smith and S. Willenbrock, Phys. Rev. D54, 6696 (1996), |arXiv:hep-ph/9604223 



T. Stelzer and S. Willenbrock, Phys. Lett. B357, 125 (1995), [arXiv:hep-ph/9505433 
G. Bordes and B. van Eijk, Nucl. Phys. B435, 23 (1995). 
S. Zhu, (2001), |arXiv:hep-ph/0109269, 



B. W. Harris, E. Laenen, L. Phaf, Z. Sullivan, and S. Weinzierl, Phys. Rev. D66, 054024 



(2002), arXiv:hep-ph/0207055 



Z. Sullivan, Phys. Rev. D70, 114012 (2004), [arXiv:hep-ph/0408049 
Q.-H. Cao, (2008). larXTv:0801.15391 

S. Frixione, E. Laenen, P. Motylinski, B. R. Webber, and C. D. White, JHEP 07, 029 (2008), 
larXiv: 0805 .30671 

J. M. Campbell, R. Frederix, F. Maltoni, and F. Tramontano, Phys. Rev. Lett. 102, 182003 
(2009), larXiv:0903.0005l 

J. M. Campbell, R. Frederix, F. Maltoni, and F. Tramontano, (2009). larXiv:0907.3933l 

J. Campbell, R. K. Ellis, and F. Tramontano, Phys. Rev. D70, 094012 (2004), 



arXiv:hep-ph/0408158 



Q.-H. Cao and C.-P. Yuan, (2004), arXiv:hep-ph/0408180 



Q.-H. Cao, R. Schwienhorst, and C. P. Yuan, Phys. Rev. D71, 054023 (2005), 



arXiv:hep-ph/0409040 



38 



[41] 

[42] 

[43] 
[44] 
[45] 

[46] 

[47] 
[48] 
[49] 
[50] 

[51] 
[52] 
[53] 
[54] 
[55] 

[56] 
[57] 
[58] 
[59] 
[60] 
[61] 
[62] 
[63] 
[64] 
[65] 



Q.-H. Cao, R. Schwienhorst, J. A. Benitez, R. Brock, and C. P. Yuan, Phys. Rev. D72, 



094027 (2005), arXiv:hep-ph/0504230 



S. Frixione, E. Laenen, P. Motylinski, and B. R. Webber, JHEP 03, 092 (2006), 



arXiv:hep-ph/05 12250 



S. Alioli, P. Nason, C. Oleari, and E. Re, JHEP 09, 111 (2009). larXiv: 0907.40761 

C. S. Li, R. J. Oakes and J. M. Yang, Phys. Rev. D55, 5780 (1997), [a7xTv:hep-ph/96TT455 



J. J. Zhang, C. S. Li, Z. Li and L. L. Yang, Phys. Rev. D75, 014020 (2007), 



arXiv:hep-ph/0610087 



M. Beccaria, C. M. Carloni Calame, G. Macorini, E. Mirabella, F. Piccinini, F. M. Renard 
and C. Verzegnassi, Phys. Rev. D77, 113018 (2008). larXTv:0802. 19941 



J. Pumplin et al, JHEP 07, 012 (2002), arXiv:hep-ph/0201195 



W. K. Tung et al, JHEP 02, 053 (2007), arXiv:hep-ph/0611254 



H. L. Lai et al, JHEP 04, 089 (2007), arXiv:hep-ph/0702268 



J. Pumplin, H. L. Lai, and W. K. Tung, Phys. Rev. D75, 054029 (2007), 



arXiy:hep-ph/0701220 



P. M. Nadolsky et al, Phys. Rev. D78, 013004 (2008). la?Xiv:0802.0007! 



Q.-H. Cao and C.-P. Yuan, Phys. Rev. Lett. 93, 042001 (2004), |aZXiv:hep-ph/040T026 
DO, V. Abazov et al, (2009). larXTv:0904.3195l 
W. T. Giele and E. W. N. Glover, Phys. Rev. D46, 1980 (1992). 

W. T. Giele, E. W. N. Glover, and D. A. Kosower, Nucl. Phys. B403, 633 (1993), 



arXiv:hep-ph/9302225 



S. Keller and E. Laenen, Phys. Rev. D59, 114004 (1999), arXiv:hep-ph/9812415 



Tevatron Electroweak Working Group, (2009). larX~iv:0903.2503l 
UA2, J. Alitti et al, Phys. Lett. B257, 232 (1991). 

ATLAS, D. O'Neil, B. Gonzalez-Pineiro, and M. Lefebvre, J. Phys. G28, 2657 (2002). 

M. Cristinziani and G. Petrucciani, Nuovo Cim. 123B, 1142 (2008), arXiv: 0808 .05651 

C. Kao, G. A. Ladinsky, and C. P. Yuan, (1993), [arXTv:hep-ph/9305270] 

C.-L. Chou and M. E. Peskin, Phys. Rev. D61, 055004 (2000). 

G. L. Kane and C.-P. Yuan, Phys. Rev. D40, 2231 (1989). 

DO, B. Abbott et al, Phys. Rev. D63, 031101 (2001), [aTxTv:hep-ex/0008024| 



C.-R. Chen, F. Larios, and C. P. Yuan, (2005), arXiv:hep-ph/0503040 



39 



[66] G. Mahlon and S. J. Parke, Phys. Rev. D53, 4886 (1996), [arXiv:hep-ph/95l2264 



[67] S. J. Parke and Y. Shadmi, Phys. Lett. B387, 199 (1996), arXiv:hep-ph/9606419 



[68] G. Mahlon and S. J. Parke, Phys. Rev. D55, 7249 (1997), arXiv:hep-ph/9611367 



10 



